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The high density of the recently discovered close-in extrasolar planet HD 149026b 

[ suggests the presence of a huge core in the planet, which challenges planet formation the- 

ory. We first derive constraints on the amount of heavy elements and hydrogen/helium 

Qh' present in the planet: We find that preferred values of the core mass are between 50 

and 8OM0, although a minimum value of the core mass is ~ 35M ffi in the extreme case 

■ of formation of the planet at > 0.5 AU, followed by late inward migration after > 1 Ga 

and negligible reheating due to tidal dissipation. We then investigate the possibility of 

. subcritical core accretion as envisioned for Uranus and Neptune. We show that a mas- 

• 1— ( , 

sive core surrounded by an envelope in hydrostatic equilibrium with the gaseous disk 
may indeed grows beyond 3OM0 provided the core accretion rate remains larger than 
~2x 1O _5 M0 yr _1 . However, we find the subcritical accretion scenario is very unlikely 
in the case of HD149026b for at least two reasons: (i) Subcritical planets are such that 
the ratio of their core mass to their total mass is above ~ 0.7, in contradiction with con- 
straints for all but the most extreme interior models of HD 149026b; (ii) High accretion 
rates and large isolation mass required for the formation of a subcritical > 35M® core 
are possible only at specific orbital distances in a disk with a surface density of dust 
equal to at least 10 times that of the minimum mass solar nebula. This value climbs 
to 30 when considering a 50M ffi core. These facts point toward two main routes for 
the formation of this planet: (i) Gas accretion that is limited by a slow viscous inflow 
of gas in an evaporating disk; (ii) A significant modification of the composition of the 
planet after gas accretion has stopped. These two routes are not mutually exclusive. 
Illustrating the second route, we show that for a wide range of impact parameters, gi- 
ant impacts lead to a loss of the gas component of the planet and thus may lead to 
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planets that are highly enriched in heavy elements. Alternatively, the planet may be 
supplied with heavy elements by planetesimals by secular perturbations. Both in the 
giant impact and the secular perturbation scenarios, we expect an outer giant planet to 
be present. Observational studies by imaging, astrometry and long term interferometry 
of this system are needed to better narrow down the ensemble of possibilities. 

Subject headings: planetary systems: formation 

1. INTRODUCTION 

The combined detection of extrasolar planets by radial velocimetry and transit photometry 
provides unique information by the measurement of both their mass and radius. Further spectro- 
scopic studies of the system in and off transit can then even inform us on chemical species in the 
planet's atmosphere, its possible evaporation, global atmospheric temperature, possible presence of 
winds, etc. 

Nine transiting extrasolar planets have been discovered so far (see www.obspm.fr/planets). 
Eight of them have radii of 1.0-1.3 Jupiter's radius (Rj), but one stands out: HD149026b has 
a radius of only 0.73 ± 0.03Rj for a mass of 0.36 ± 0.03 Jupiter's mass (= 110 ± 10M e ) and 
an orbital distance of 0.042 AU (Sato et al. 2005; Charbonneau et al. 2006). The parent star is 
believed to be 1.3 M Q , and is metal-rich with [Fe/H] = 0.36±0.05. Contrary to the other transiting 
extrasolar planets, believed to be mainly formed with hydrogen and helium (Burrows et al. 2000; 
Bodenheimer et al. 2001; Guillot & Showman 2002; Laughlin et al. 2005; Guillot 2005; Baraffe et 
al. 2005), HD149026b is clearly made of a significant fraction of heavy elements. 

Indeed, based on evolution models, Sato et al. (2005) and then Fortney et al. (2006) have shown 
the planet has a huge rock/ice core of ~ 70M e . As noted by Sato et al. (2005), such a big core points 
toward formation of the planet in the core accretion scenario (e.g., Mizuno 1980; Bodenheimer & 
Pollack 1986), rather than in the disk instability scenario (e.g., Boss 1997). However, HD149026b 
may even challenge the core accretion scenario, because the core mass is far larger than that of 
Jupiter or Saturn, and the "canonical" critical core mass as derived by Mizuno (1980) is around 
10M e . 

The purpose of this paper is to constrain the structure of HD149026b, to examine how the 
planet may have formed based on these constraints, and, hopefully, to stimulate further studies 
of the formation of close-in giant planets. In section 2, we attempt to constrain the structure 
of HD 149026b based on various evolution models. In section 3, we then examine whether these 
constraints are compatible with the properties of growing protoplanetary cores embedded in a 
gaseous circumstellar disk. On these bases, various possibilities for the formation of this planet 
are explored in section 4. Section 5 summarizes the results and proposes two possible formation 
scenarios that account for the properties of HD149026b. 
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2. PRESENT STRUCTURE: HOW MUCH HEAVY ELEMENTS? 

As already evident from the measurements and from previous models of the evolution of the 
planet (Sato et al. 2005; Fortney et al. 2006), HD149026b is surprisingly small, and should possess 
either a big core or a large amount of heavy elements in its interior. In this work, we attempt to 
constrain the core mass by using a variety of evolutionary models beyond those explored by Sato 
et al. and Fortney et al. In doing so, we will show that HD149026b indeed contains a significant 
amount of heavy elements, in apparent contradiction with current models of the formation of giant 
planets. 

2.1. Standard Models 

The evolution of a giant planet is essentially governed by a Kelvin-Helmholtz cooling and 
contraction similar to a stellar pre-main sequence evolution but slightly modified by degeneracy 
effects, and by the intense stellar irradiation that slows the cooling through the growth of a deep 
radiative region that can extend to kbar levels (Guillot et al. 1996; Guillot 2005). 

The level of the irradiation received by the planet is calculated from the observations by Sato et 
al. (2005): We adopt values of the stellar radius (i?* = 1.45 R ) and temperature (T* = 6150 K) and 
the orbital distance of 0.046 AU to derive zero-albedo equilibrium temperature T c * q = 1740 K. We 
use this value as a minimum for the atmospheric temperature at the 1 bar level (the corresponding 
temperature found by Fortney et al. (2006) is between 2000 and 2200 K, and using the model of 
Iro, Bezard, h Guillot (2005), we find a value of 1980 K). We use the 1 bar level as outer boundary 
condition for the evolution models (see Guillot 2005). 

Standard models are calculated on the basis of structure consisting of a central rock/ice core 
and of an envelope of solar composition, and further assuming that the planet has either formed 
in situ or moved rapidly to its present location, so that throughout its evolution it has received a 
constant stellar heat flux to slow its cooling. 

For the core, we use the equations of state for "rocks" and "ices" obtained by Hubbard & 
Mar ley (1989). For the envelope, we use the H-He EOS from Saumon, Chabrier, & Van Horn 
(1995). The structure of the core with Hubbard & Marley (1989) 's EOS is independent of its 
temperature, but this simplification is reasonable because it has much smaller effects than changes 
in the (unknown) core composition. Although it has little effect on the models, we assume that the 
core luminosity is due both to radioactive decay (with / c hondritic ~ 10 20 ergs -1 M^ 1 ) and cooling 
(assuming that in the core the temperature is uniform and the specific heat is c v « 10 7 ergg -1 K -1 ). 
Internal opacities are calculated from Rosseland opacity tables provided by Allard et al. (2001). 

The calculations are made for various core masses and core compositions and compared to the 
observational constraints in Fig. 1. As discussed by Sato et al. (2005) and by Fortney et al. (2006), 
the small planetary radius implies that it contains a large fraction of its mass in heavy elements. 
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Difference in composition of the core — whether these heavy elements are mostly "ices" (a mixture of 
water, ammonia and methane in their [unknown] high pressure, relatively high temperature form) 
or "rocks" (a mixture of refractory materials including mostly silicates) — affects the radius of the 
planet by 10 to 20%. For comparison, the interiors of Uranus and Neptune are consistent with 
being mostly made of "ices" (e.g., Podolak, Podolak, & Mar ley 2000). Figure 1 thus confirms the 
need for a substantial amount of heavy elements in HD149026b when adopting a standard evolution 
scenario. Quantitatively, this implies core masses between 41 and 83 in good agreement with 
Sato et al. (2005) and Fortney et al. (2006). In particular, the latter yields values between 60 and 
93 M ffi . The difference is due to the fact that in order to constrain more strictly the minimum mass 
of heavy elements present in the planet, we used a low value of the atmospheric temperature. 



2.2. The Cold Storage Hypothesis 

We now investigate the hypothesis that the planet may have cooled far from its parent star 
before suddenly being sent into the presently observed 0.046 AU orbit. The reason for this sudden 
event is to be determined, but could be due to dynamical planet-planet interactions (see section 
4.1.2). In terms of evolution, the planet is allowed to cool more rapidly during the time-period 
when it is far from the parent star. The sudden inward migration indeed yield an expansion of the 
planet but it may be limited to the outer layers and thus be relatively small (Burrows et al. 2000). 
We want to investigate whether the present mass and radius are compatible with a relatively small 
core. 

We thus calculate evolution models similar to those in the previous section, but using T eq = 
100 K and an atmospheric model based on simplified radiative transfer calculations of an isolated 
atmosphere (see Saumon et al. 1996; Guillot 2005). Based on these calculations, we obtain generally 
a much faster contraction of the planet; a less than 30 M® core can account for the observed 
radius. However, we have to include reheating of the outer shell of the envelope. This shell 
is approximatively defined by the region whose temperature is lower than the new equilibrium 
temperature, 1740 K. Based on our calculations, this corresponds to a pressure level of ~10kbar. 

A simple estimation of the reheating timescale of the outer shell suggests that we have to 
include this effect. Using equations of radiative diffusion and energy conservation, we estimate that 
the reheating timescale of a layer of pressure P, temperature T, opacity k and heat capacity c p is 
to first order independent of its initial temperature and equal to: 

^ 3c p kP 2 . . 

Trad ~4^?T^ (lj 

where a is the radiation density constant, c the speed of light, and g the planet's gravity (assumed 
uniform). This expression is consistent with direct simulations by Iro, Bezard, & Guillot (2005). 
This timescale is about 10 8 years for k rs lcm 2 g _1 , P rs lOkbar, g k, 1550 cm s -2 , T k, 2000 K and 
Cp R^ 8 X 10 7 ergg 1 K 1 , which implies that the reheating is relatively fast and should involve most 
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of the shell that we defined. 

However, the expansion due to the reheating is relatively small. The P 2 dependence ensures 
that most of this outer shell will be affected, and it will expand by a radius estimated to be: 

AR ss H p AlnP, (2) 

where H p « TZT/ fig is the new pressure scale height and A In P measures the depth of the expanding 
shells in pressure units. Based on our simulations, a rough estimate is that the planet should expand 
by ~ 4 x 10 8 cm (= 0.06-Rj) when considering levels between lOkbar and 1 bar. 

Note that this estimate does not include factors such as global reheating due to tidal dissipation 
caused by the necessary circularization of the planet's spin and eccentricity. If such dissipation 
occurs sufficiently deep in the planet, this will lead to an expansion of the planet's deeper layers 
that could completely erase the planet's cold evolution phase (see Guillot 2005, and references 
therein). Thus, this very possibility that a planet can contract more efficiently far away from the 
star and then be brought in with little consequence on its radius is very much in question. We 
choose however to consider it because it gives a sense of the robustness of the conclusion that 
HD 149026b contains a large fraction of heavy elements. 

Table 1: Derived constraints on HD149026b's composition. MnHe and Mz are masses of H/He and 
heavy elements, and M p = Mhhc + Mz is total planetary mass. 





M z (M ffi ) 


Mhhc 




M z /M p 


In situ evolution 








ices 


50 - 83 


27- 


54 


0.52-0.72 


rocks 


41 - 69 


38- 


64 


0.43-0.60 


Cold storage scenario 






ices 


37-77 


34- 


66 


0.39-0.67 


rocks 


33-61 


46- 


76 


0.33-0.53 


Total 












33-83 


27- 


76 


0.33-0.72 



Adding AR (~ 4 x 10 8 cm) to the radius obtained by the evolution calculation of the isolated 
planet (not including the global reheating) , we derive masses of heavy elements that reproduce the 
observed radius and show these in Table 1. Even in the case of a long storage (> 10 8 years) of 
the planet in a cold environment (corresponding to more than 10 AU), a large amount of heavy 
elements (strictly more than 33 M®) is required to reproduce the observed radius of HD149026b. 
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2.3. Heavy Elements: In the Core or in the Envelope? 

We have thus far assumed all heavy elements to lie within a central core. In order to see how 
the results are dependent on this assumption, we compare our previous model with a 60M ffi ice 
core to a similar model with a 30 M e ice core and an envelope which is enriched by 30 M® of ices. 
For the ices in the envelope, we use the EOS described by Saumon & Guillot (2004). 

Figure 2 shows the results of the calculations in the two cases: (i) with unchanged opacities 
(i.e., as calculated for a solar-composition mixture) and (ii) with 30 times larger opacities to mimic 
the effect of the enrichment of heavy elements in the envelope. Basically, when using an unchanged 
opacity table (the dotted line), there is very little difference between a planet that has all its heavy 
elements in the core and a planet that has them mixed throughout its envelope. However, large 
differences arise when the opacities are affected proportionally to the amount of heavy elements 
that are mixed (the dashed line). In that case, the cooling and contraction timescale, which is 
dominated by radiative transport in the outer radiative zone, becomes long, and prevents the rapid 
contraction of the planet. 

These results confirm again that our estimates of the amount of heavy elements in HD149026b 
are lower limits and that the planet indeed must contain a significant amount of heavy elements. 
If most of the heavy elements are located below the radiative zone, which is not unlikely, then we 
expect our constraints to be relatively accurate. On the other hand, if the material is mixed in the 
outer radiative zone, we expect that its interior contains more heavy elements by up to ~ 20M ffi 
than calculated here. 

2.4. Impact of the Various Parameters 

We have used a rather simplified approach for the calculations of the interior structure. On the 
other hand, reality is without doubt more complex. However, this should not affect significantly 
the global constraints that are derived in Table 1. First, the two extreme compositions used for 
the heavy elements ("ices" and "rocks") ensure that most variations due to improved EOS, for 
example, are likely to fall in between the range of values that are considered. 

Our external boundary condition is extremely simplified. Although it agrees with more detailed 
atmospheric models, one has to account for the fact that the atmospheric temperatures may be 
higher, in particular if the atmosphere is itself enriched in heavy elements (see Fortney et al. 2006). 
A higher atmospheric temperature leads to a slightly larger radius, everything else being the same. 
However, given the decrease of the opacity with increasing temperature in that regime, the radius 
increase is very limited. 

Other factors point toward a slightly larger amount of heavy elements than calculated here: 
(i) as discussed, the likely increase of the opacities in an enriched envelope; (ii) the presence of any 
other energy source such as the one that is required to reproduce the radius of HD209458b; (iii) 
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thermally-dependent EOS for the core. Again, given that the core mass is already quite large, we 
expect these factors not to change the masses of heavy elements that are inferred above. 

3. A SUBCRITICAL CORE ACCRETION SCENARIO? 

The structure of HD 149026b as derived in the previous section is certainly puzzling: It contains 
as much hydrogen and helium as Saturn, but at least twice as much heavy elements. In other words, 
the ratio of heavy elements to hydrogen and helium is intermediate between those of Saturn and 
Uranus/Neptune. In the core accretion scenario, when core mass exceeds a critical value (called 
the critical core mass), the Kelvin-Helmholtz contraction of the envelope takes place, resulting in 
substantial disk-gas accretion onto the planet (e.g., Mizuno 1980; Bodenheimer & Pollack 1986). 
A widespread idea is that Jupiter and Saturn experienced the substantial gas accretion beyond 
the critical core mass, while the cores of Uranus and Neptune remained subcritical because of 
progressive accretion of planetesimals until the gaseous circumstellar disk disappeared (Pollack et 
al. 1996). In this section, we investigate whether HD149026b may have formed in a way similar to 
what envisioned for Uranus and Neptune. 

3.1. Critical Core Mass and Gas Accretion Rate 

Figure 3 shows the critical core mass (M CiCr jt) as a function of core accretion rate for several 
different choices of distance from the parent star and local density of disk gas. In the numerical 
simulations we have used the H-He EOS from Saumon et al. (1995), the grain opacity from Pollack, 
McKay, & Christofferson (1985), and the gas opacity from Alexander & Ferguson (1994). The 
integration method is basically the same as that of Ikoma et al. (2000) who used simpler forms of 
EOS and gas and grain opacities. 

The input physics used here is somewhat different from that used in section 2. First we include 
grain opacity that was not included in the simulations in section 2, because an accreting envelope 
contains low temperature parts in which grain opacity is dominant, unlike the fully-formed planet 
(T > 1700 K). Although the fully-formed planet also contains such low temperature parts in the 
cold storage model, the planet is almost fully convective. Second the gas opacity used in this 
simulation is different from that used in the previous section. This difference has little influence on 
the results obtained in this section, because most of the part of the envelope where gas opacity is 
dominant relative to grain opacity is convective. Finally, although constant core density is assumed 
here, the critical core mass is known to be insensitive to core density (Mizuno, Nakazawa, & Hayashi 
1978). 

Figure 3 shows a large subcritical core can be formed if core accretion rate is high. High core 
accretion rates stabilize the envelope, yielding progressively larger critical core masses (Ikoma et al. 
2000). This, however, remains true up to a point for which the envelope becomes fully convective; 



-8- 



the structure of a fully convective envelope is almost adiabatic and thus independent of energy flux 
supplied by incoming planetesimals (i.e., core accretion rate). At that point, the critical core mass 
reaches its maximum value. In this case, this maximal critical core mass depends on the local disk 
conditions (i.e., the local density and temperature of disk gas and the distance from the parent 
star) (Wuchterl 1993; Ikoma, Emori, & Nakazawa 2001). In the minimum-mass solar nebula (MSN) 
proposed by Hayashi (1981), the entropy of the disk gas decreases as the distance to the parent 
star decreases. That is why the critical core mass is smaller at locations closer to the parent star: 
In a denser gas disk, the critical core mass is smaller for the same reason (Ikoma et al. 2001). As a 
result, formation of a large subcritical core of 50-80 M ffi requires not only high core accretion rate, 
but also a location not too close to its parent star and not too massive gaseous disk. 

Although a large core of 50-80 can, in principle, be formed by subcritical growth as shown 
above, the ratio of the critical core mass to the corresponding total (core + envelope) mass suggests 
that the core mass of HD 149026b was likely to be supercritical. Figure 4 shows the ratio of core 
mass to planetary total mass as a function of core mass (normalized by the critical core mass) for 
three different sets of core accretion rate and distance from the parent star. This ratio is found to 
be rather insensitive to values of the parameters; it is always ~ 0.7 at the critical core mass. The 
evolution model of HD 149026b obtained in section 2 shows that the current ratio of core mass to 
planetary mass is smaller in most cases (see Table 1). Although there are a few cases where the 
current ratio is as large as ~ 0.7, the probability of the formation of such planets seems to be low. 
From Fig. 4 we also find the duration during which the ratio of core mass to planetary mass is 
around ~ 0.7 is quite short in the total core formation time. It follows from this consideration that 
HD 149026b is likely to have experienced disk-gas accretion due to the Kelvin-Helmholtz contraction 
of the envelope. 

The rate of the gas accretion determined by the Kelvin-Helmholtz contraction of the envelope 
beyond the critical core mass — the gas accretion is hereafter called the KH gas accretion — depends 
strongly on the critical core mass. Figure 5 shows the typical timescale (r g ) for the KH gas accretion 
defined by 

,crit ,crit / \ 

r s = c 5 t ' ( 3 ) 



where G is the gravitational constant, M CjCr i t is the critical core mass, M ejCr i t is the corresponding 
envelope mass, R com is the outer radius of the inner convective region, and L is the luminosity 
(Ikoma et al. 2000). The numerical factor, C, is chosen to be 3/2 based on the numerical simulations 
by Ikoma & Genda (2006). In cases where core accretion is suddenly stopped (corresponding to 
core isolation; see below), C should be 1/3 (Ikoma & Genda 2006). On this timescale the envelope 
mass increases by a factor of e. The timescale for the KH gas accretion is approximately expressed 

by 

r g ~5xl0 10 (^)" 3 ' 5 years (4) 
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in this range . This relation implies that if the core mass is small, the timescale of the KH gas 
accretion is much longer than the disk lifetime (~ 10 6 -10 7 years). 

For discussion in section 3.2, we have also drawn lines representing the time for formation of 
a critical core with constant core accretion rate, 

- Mi i a* 

7"c,acc,cnt — ^ \ M c =M CjCrit > {?) 

for three different sets of the parameters. 



3.2. Subcritical Core Accretion 

We examine the subcritical core accretion to know how large core can be formed before the 
onset of the KH gas accretion. Core mass is limited by the isolation mass that is determined by 
the total mass of solid components in the core's feeding zone, unless planetesimals are supplied 
from outside the feeding zone by migration of planetesimals or the core itself. Merging between 
the isolated cores is not likely until disk gas is severely depleted (Iwasaki et al. 2002; Kominami 
& Ida 2002), because the eccentricity damping due to disk-planet interaction (Artymowicz 1993; 
Ward 1993) is too strong to allow orbital crossing of the isolated cores. The core isolation mass is 
thus given by (Ida & Lin 2004a) 

where a is semimajor axis of the core. The scaling factor / g and /d express enhancement of disk 
surface density of gas (£) and dust (a) components relative to those of MSN, defined by 

E " 2400 ^(ilu)" 3/2g/cm2 m 

( a \-3/2 2 \') 
a = 10/d??icc t~ttt g/cm . 



.1AU, 

In the solar metallicity cases, the values of these factors may have a range between 0.1 and 10 (Ida 
& Lin 2004a). For metal-rich disks, the range of /d may shift to larger values by a factor of 10' Fe / H l, 
while that of / g does not change (Ida & Lin 2004b) . 7?i cc expresses an enhancement of disk surface 
density due to ice condensation, which is 3-4 beyond an ice line at ~ 3 AU. Equation (6) shows 
large cores form at large a, in particular beyond the ice line, in disks with large f^. As seen from 
Fig. 3, smaller core accretion rate yields smaller critical core mass, which implies core isolation 
triggers the KH gas accretion. For a core of M c not to be isolated, that is, M c < M Cj i so , 

/d>30 — H ^ice ( TTTF ) IT 1 • (8) 



30M ffi y ,1CC UALV V M ©/ 



1 The form of equation (4) is different from that of equation (19) of fkoma et al. (2000). The difference comes 
mainly from the difference in grain opacity (see Ikoma & Genda 2006, for the detailed discussion) . 
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Next, we consider the condition for a core to avoid the KH gas accretion until its isolation. 
Figure 3 shows that M CjCr i t is expressed approximately by 

M ccrit ~ 10 [ 7 Mc - | M ffi , (9) 

^10- 7 M e /yr J ®' V ' 

except for fully convective cases. Note that consideration of fully convective cases does not change 
our conclusion (see below). The condition for a core to avoid the KH gas accretion is M c < M CjCr it- 
Using T Ciacc = M c /M c and Eq. (9), this condition is rewritten by (see Fig. 5) 

(. M " 



T c ,acc < T CiacCiCri t ~ 1 x 10 6 J years. (10) 

For cores of 30, 45 and 70M®, T Ciacc must be shorter than 1 x 10 6 years, 2 x 10 5 years, and 3 x 10 
years, respectively. The core growth timescale due to planetesimal accretion is given by (Ida k, Lin 
2004a) 

r ~ 9 3 x 10 5 V 1 (—) 27/10 (^~) (^A (— "1 2/15 years (11) 

r c , acc _ 9.3 x 10 / d / g ^ e ^ 1AU J V30M e y \M Q J \l^g) Y ^ {U) 

where m is the mass of a planetesimal. The dependence on f g appears because damping of velocity 
dispersion of planetesimals due to aerodynamical drag is taken into account. Substituting Eq. (11) 
into Eq. (10), we have 

This inequality indicates that large /d is needed to keep core accretion rate high enough to avoid 
the KH gas accretion. 

For subcritical core accretion, both Eqs. (8) and (12) must be satisfied. For M* = 1.3Mq and 
[Fe/H]=0.36 (/d// g = 2.3), which correspond to HD149026, the conditions are plotted in Fig. 6. 
Subcritical core accretion could be possible far from the parent star in relatively massive disks for 
30M© cores. However, for 45 and 70M e cores, extremely massive or extremely metal-rich disks 
with /d > 30 — 50 are required. Even for those large values of /d, since the critical core mass reaches 
its maximal value in cases of high core accretion rate, subcritical formation of 45 and 70 M e is 
impossible in some cases (see Fig. 3). Although formation of a 3OM0 core may not be very difficult 
in relatively metal-rich disks, 30Me is the lowest estimated value for Mz of HD 149026b in the 
cold storage model; 50-80Me is most likely. Therefore, heavy elements must be supplied after 
substantial accretion of the envelope. This issue is discussed in section 4.1. 



4. SEVERAL POSSIBILITIES FOR THE FORMATION OF HD149026b 

The investigations in sections 2 and 3 have revealed that we have to explain at least three 
properties of HD149026b: its small orbital radius, its high metalicity, and its moderate amount of 
hydrogen and helium. 
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A possible scenario that explains the proximity of the planet to its parent star is the in situ 
formation or migration of the planet. The in situ formation is however unlikely, because the core 
isolation mass is so small at ~ 0.04 AU (see eq. [6]) that the planet can get only a tiny amount of 
hydrogen/helium within the disk lifetime because of too long timescale for the KH gas accretion 
(see eq. [4]). 

Section 3.2 has shown cores of > 301M® can form only far from the parent star in reasonably 
massive and/or metal-rich disks. Except for the cold storage model (section 2.2), residual heavy 
elements of 20-50M® must be supplied to the planet. In section 4.1, we examine two possible 
scenarios; supply of heavy elements (planetesimals or another planet) during gas accretion (§ 4.1.1) 
and after disk gas depletion (§ 4.1.2). To account for the moderate amount of hydrogen and helium 
of HD149026b, we examine limited supply of disk gas (§ 4.2.1) and loss of the envelope gas (§ 4.2.2). 

4.1. Supply of Heavy Elements 

4-1.1. During gas accretion 

First we briefly discuss the possibility that the large amount of heavy elements in HD149026b 
were supplied during substantial accretion of the envelope that happened far from the parent star. 
An increase in the planet's mass due to gas accretion expands its feeding zone. If the planet 
can efficiently capture the planetesimals in the feeding zone, supply of heavy elements during gas 
accretion is possible (Pollack et al. 1996; Alibert, Mordasini, & Benz 2004). However, the cross 
section of gravitational scattering is much larger than that of collision in outer regions. Most 
planetesimals are thus scattered before accretion onto the planet. The gravitational scattering 
combined with eccentricity damping by disk-gas drag ends up opening a gap in the planetesimal 
disk: This phenomenon is called "shepherding" (e.g., Tanaka & Ida 1997) 2 . Although fast migration 
or growth of the planet can avoid the gap opening, unrealistically fast migration or growth may be 
required for a Saturn- mass planet to avoid the gap opening (Tanaka & Ida 1999; Ida et al. 2000). 
Thus, supply of a large amount of planetesimals may be unlikely to happen during gas accretion. 

However, much more detailed studies are needed for this possibility. If the shepherding com- 
pletely prevented the planetesimal accretion throughout the gas accretion phase, the well-known 
enrichment of heavy elements in the envelopes of Jupiter and Saturn would be unable to be ac- 
counted for. A possibility is accretion of small fragments. Small fragments are produced by dis- 
ruptive collisions of planetesimals in the vicinity of a large core (Inaba, Wetherill, & Ikoma 2003). 
Since atmospheric-gas drag is strong for the fragments, they are accreted onto the core (Inaba & 
Ikoma 2003). Because these studies are limited to the subcritical core case, we need to extend 
the calculation to the case of supercritical cores, in which the planet mass rapidly increases and 



2 The shepherding by isolated mean-motion resonances is called "resonant trapping" (e.g., Weidenschilling & Davis 
1985; Kary, Lissauer, & Greenzweig 1993). 
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non-uniformity of disk gas is pronounced. 

4-1-2. After disk gas depletion 

We next explore the possibility of enrichment of heavy elements in HD149026b (i.e., bombard- 
ment of another planet or planetesimals) after substantial accretion of the envelope. 

If the late bombardment happened, it is likely to have occurred after the planet migrated to 
its current location: The close-in planet collided with another planet (s) or planetesimals that had 
been gravitationally scattered by an outer giant planet(s). At ~ 0.05 AU the planet's Hill radius 
is only several times as large as its physical radius, while the former is much larger (by 100-1000) 
than the latter in outer regions (> 1 AU). That means the ratio of collision to scattering cross 
sections is much larger at ~ 0.05 AU relative to in outer regions. The scattering cross section 
is further reduced by high speed encounter. The impact velocities (t>i mp ) of the scattered bodies 
in nearly parabolic orbits are as large as the local Keplerian velocity at ~ 0.05 AU, which is a 
few times larger than surface escape velocity of the inner Saturn- mass planet (f CS c)- The collision 
cross section is, in general, larger than the scattering cross section for v\ mp > v csc . Note that the 
shepherding does not work for highly eccentric orbits. 

To know how efficiently such scattered bodies collide with the close-in giant planet, we perform 
the following numerical simulation. We consider an inner planet of 0.5Mj in a circular orbit at 
0.05 AU and a planetesimal (test particle) or another giant planet of 0.5Mj in a nearly parabolic 
orbit of e ~ 1 and i = 0.01 with initial semimajor axis of a = 1 AU. For each initial pericenter 
distance (q), 100 cases with random angular distributions are numerically integrated by 4th order 
Hermite integrator for 100 Keplerian periods at 1 AU. Then, collision probability with the inner 
planet (-P C oi)> that with the parent star (-P c * ol ), and ejection probability (-P c jc) are counted. The 
residual fraction, 1 — (P co i + P c * ol + _P e jc), corresponds to the cases in which the planetesimal/planet 
is still orbiting after 100 Keplerian periods. We also did longer calculations and found that P co \/P e j c 
is similar, although individual absolute values of P co i and P e jc increases. The results are insensitive 
to a as long as o > 0.05 AU, because velocity and specific angular momentum of the incoming 
planetesimal/planet are given by v ~ y/ 2GM* / q and L ~ \J2GM*q that are independent of initial 
semimajor axis a of incoming bodies. 

Figure 7 shows the probabilities of the three outcomes of encounters between the inner planet 
and a planetesimal (panel a) or another giant planet (panel b). As shown in Fig. 7, when 0.01AU < 
q < 0.06AU, the incoming bodies closely approach the inner planet. Although some fraction of 
them are ejected, the comparable fraction results in collision with the inner planet. In particular, 
P co i > Pcjc in the planet-planet case (panel b), in which the ratio of geometrical cross section to 
Hill radius is larger than in the planet-planetesimal case (panel a). If q is smaller than or close to 
the parent star's physical radius (which is 0.01 AU in the calculations here), most of the incoming 
bodies hit the parent star rather than collide with the inner planet. If q > 0.06 AU, they do not 
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closely encounter with the inner planet, so that both collision and ejection probabilities almost 
vanish (P co i,-P e jc < 1)- 

The results shown in Fig. 7 demonstrate that efficient supply of heavy elements to the inner 
giant planet requires a very limited range (close to unity) of eccentricities of the incoming bodies; 
for example, 0.94 < e < 0.99 for a = 1 AU, since q = a(l — e). In the case of chaotic scattering 
by outer giant planets, the probability to acquire q in such a narrow range would be very small 
(e.g., Ford, Havlickova, & Rasio 2001). On the other hand, if the scattering comes from secular 
perturbation, e is increased secularly while a is kept constant. Accordingly, q decreases secularly, 
and the incoming bodies eventually collide with the inner planet. 

For example, the following scenario is possible. More than three giant planets in addition to the 
inner one are formed in nearly circular orbits. They start orbit crossing on timescales longer than 
their formation timescales (Rasio & Ford 1996; Weidenschilling & Marzari 1996; Lin & Ida 1997; 
Marzari & Weidenschilling 2002). In many cases, one giant planet is ejected and the residual giant 
planets remain in stable eccentric orbits after the orbit crossing (Marzari & Weidenschilling 2002). 
In about 1/3 cases, q of a giant planet decreases secularly to < 0.05 AU due to the Kozai effect 
(Kozai 1962) during the orbit crossing (Nagasawa et al. 2006). This scenario requires formation of 
at least four giant planets. 



4.2. Origin of the Small Amount of Hydrogen and Helium 

4-2.1. Limited disk gas supply 

Gas accretion onto a core can be truncated by gap opening or global depletion of disk gas. A 
gap may be formed, if both the viscous and thermal conditions are satisfied (e.g., Lin & Papaloizou 
1985, 1993; Crida, Morbidelli, & Masset 2006). The former is that the torque exerted by the planet 
overwhelms viscous torque, and given by (Ida & Lin 2004a), 

^Hiluf (SH (13) 

where M p is the planet mass and a is the parameter for the a-prescription for disk viscosity 
(Shakura & Sunyaev 1973). The latter condition is that the Hill radius becomes larger than disk 
scale height, and given by 3 (Ida & Lin 2004a) 

/ a \ 3/4 / M* \ 

M p >M p , lt = 400(— ) (— }Me. (14) 

assuming an optically thin disk in which T = 280(a/lAU)~ 1 / 2 K. Hence, the actual truncation 
condition may be the thermal condition and it is very unlikely that gas accretion is truncated by 



3 Here we assume the Hill radius (rH,c) is equal to 1.5 times disk scale height (h), roughly taking into account late 
gas accretion with a reduced rate after rH, c = h, following Ida & Lin (2004a). 
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the gap opening at M p ~ 110M e (which corresponds to HD149026b) far from the parent star, 
as long as a sufficient amount of disk gas remains. Furthermore, since the gap is replenished by 
viscous diffusion, gas accretion may not completely be truncated by the gap opening (D'Angelo, 
Kley, & Henning 2003). 

If a sufficient amount of disk gas remains, one way to truncate gas accretion at relatively 
small planetary mass is that the planet migrates to the vicinity of its parent star before the planet 
accretes a large amount of gas. Both M PiV ; s and M p t h are small in the inner regions (Eqs. [13] and 
[14]). However, this is unlikely, even if type II migration occurs when M p ~ 30M®. The timescale 
for the KH gas accretion < 3 x 10 4 years for a core of > 30M ffi (Fig. 5), while the timescale for the 
migration is much longer (~ 10 6 -10 7 years). Furthermore, it is not clear if the gap opening can 
stop gas accretion completely. 

If disk gas is globally depleted when M c reaches M CiCr it, gas accretion onto the core can be 
limited by viscous diffusion of disk gas, not by the Kelvin-Helmholtz contraction of the envelope 
(e.g., Guillot h Hueso 2006). That could be possible to account for the small amount of H/He. 
However, whether such a small amount of disk gas can bring the planet to the vicinity of the parent 
star should be examined. If it does not work, different migration mechanism such as gravitational 
scattering during orbital crossing of giant planets is required. 

4-2.2. Loss of envelope gas 

Another way to explain the small amount of the H/He gas of HD 149026b is loss of the envelope 
gas. There are three possibilities for loss of the envelope gas; photoevaporation driven by incident 
UV flux from the parent star, the Roche lobe overflow, and impact erosion by a collision with 
another gas giant planet. 

Photoevaporation process for gas-rich planets is normally faster for smaller planetary masses, 
because less massive planets are more expanded and their envelope gas is more loosely bound 
(Baraffe et al. 2005). This means the envelope quickly disappears once the evaporation occurs, 
so that it should be relatively rare to observe a planet at a stage when a relatively small amount 
(30-60 M e ) of envelope gas remains. However, this possibility is not excluded at present, because 
Z-rich planets such as HD149026b can be more compact and their envelope gas is not necessarily 
more loosely bound (e.g., Baraffe et al. 2006). 

Envelope gas can also be lost by the Roche lobe overflow. When the Roche lobe is filled with 
the envelope gas, the gas overflows through the inner Lagrangian point to the inner disk. The Roche 
lobe overflow takes place if a planet is very close to its parent star because of inflation by stellar 
tidal heating as well as reduction of the Roche lobe (e.g., Trilling et al. 1998), and subsequent tidal 
inflation instability of the envelope (e.g., Gu, Lin, & Bodenheimer 2003). The spilled gas from the 
Roche lobe first goes to the inner disk and is eventually going to fall on the parent star by the 
gravitational interaction with the planet. At the same time, the planet gains angular momentum 
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as the counteraction and migrates outward. Because of the increase in the orbital radius, the 
evaporation process slows down. Thus a state of a 30-60 M ffl envelope could be possible. This 
possibility is, however, in question because it is unclear if sufficient angular momentum is lost. 

Sato et al. (2005) suggested that collision between two giant planets could account for the 
internal structure of HD149026b. If cores with individual mass of ~ 30M® are merged while 
significant amounts of their envelopes evaporate, a planet similar to HD149026b is formed. Here 
we examine if significant loss of the envelope gas occurs through SPH simulations of a collision 
between two giant planets. 

High impact velocity is expected. Suppose that one giant planet orbits at 0.03-0.05 AU and 
another giant planet in a highly eccentric orbit is sent from an outer region (for the mechanism to 
send planets, see section 4.1.2). Since the impactor has high orbital eccentricity (e ~ 1), collision 
velocity is similar to local Keplerian velocity at 0.03-0.05 AU around a 1.3M Q star, 150-190 km/s. 
This velocity is 2-3 times as high as the two-body surface escape velocity. Since specific angular 
momentum of the impactor is given by L = \J GM*a(l — e 2 ) ~ y/2GM*q, the impactor and the 
target have similar L. Hence, the semimajor axis of the merged planet would be similar to that of 
the target planet unless envelope evaporation changes specific angular momentum significantly. 

Previous SPH simulations of a collision between Mars-size or Earth-size rocky planets (e.g., 
Canup 2004; Agnor & Asphaug 2004) show almost no loss of material occurs by a low-velocity 
impact (vimp) of 1.0-1.5 times as high as the two-body surface escape velocity (t> esc ) defined by 
fg SC = 2G(M t +Mi)/(R t + Ri), where M and R are the planetary mass and radius and subscriptions 
t and i represent the target and the impactor. Collisions with higher impact velocity lose larger 
amounts of material; for example, 40% of the total mass is lost when v imp ~ 2.5f csc in their 
simulations. For gas giant planets, we will have qualitatively similar results, but can expect more 
significant loss of envelope gas because the envelope is less tightly bound compared to rocky planets. 

In order to simulate the collision between two gas giant planets, we use a SPH method with the 
modified spline kernel function (Monaghan Sz Lattanzio 1985) and artificial viscosity (Monaghan & 
Gingold 1983). Our SPH simulation code was checked to reproduce the results of Agnor k, Asphaug 
(2004) for Mars-size planets using the Tillotson EOS. In this paper, we assume that the gas giant 
planets are composed only of ideal gas with the adiabatic exponent of 2; this value reproduces 
the present structure of Jupiter. We consider two hydrostatically equilibrated planets (each being 
composed of about 30,000 SPH particles) with Jupiter's mass and radius for initial conditions. We 
perform the simulations of head-on collisions with various impact velocities (1.0, 1.5, 2.0, 2.25, 2.5, 
and 3.0v csc )- In our calculations, v csc is 60 km/s. We calculate the gravitational force for each SPH 
particle using the special-purpose computer for gravitational N-body systems, GRAPE-6. 

Figure 8 shows the mass fraction lost by collisions with various impact velocities. Higher 
impact velocity results in higher loss fraction. Almost no escape occurs for a low-velocity impact 
(vimp ~ v csc ). When the impact velocity is higher than 2.5t> eS c, almost all SPH particles (~ 90%) 
are lost. 
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Although we are unable to know exactly the loss fraction of the core with our current simu- 
lation, we can estimated it approximately by tracking the motion of the SPH particles which arc 
initially located in inner part of the planet. Figure 8 also shows the loss fractions of the inner parts 
that are initially located inside 1/2 and 1/4 of the planetary radius. The inner part of the planet 
tends to remain in the merged planet. For example, in the case of 2.0v esc , about 50% of materials 
is totally lost, but almost no escape occurs for the material inside of 1/4 initial radius. A head-on 
collision of two gas giant planets with v nnp ~ 2-2.5t> csc possibly results in the considerable depletion 
of envelope and merging of core. 

The strongest point of this model for envelope loss is that we need no additional process 
to supply heavy elements to the planet. The collision increases Mz up to ~ 60M® as well as 
decreases MnHe! as we showed in section 3.2, it is possible to form cores of ~ 3OM0 by ordinary 
planetesimal accretion in relatively massive disks, in particular, metal-rich disks. However, the 
numerical simulations performed here is still preliminary. Since this mechanism may be promising, 
we should perform detailed simulations in the future. 

4.3. Other Possibilities 

Another possibility of formation of a metal-rich giant planet is the formation in a disk with 
originally high dust/gas ratio. In general, dust grains migrate inward in a gas disk to change 
the local dust/gas ratio (e.g., Stepinski & Valageas 1997). If the dust/gas ratio is enhanced in 
planet-forming regions at < 10 AU, the metal-enriched envelope could be formed without any 
planetesimal/planet supply. The high dust/gas ratio also makes possible planetesimal formation 
from dust through self-gravitational instability against Kelvin-Helmholtz instability (Youdin & Shu 
2002) and avoids type-I migration of terrestrial planets (Kominami, Tanaka & Ida 2005). However, 
it is not clear if the dust /gas ratio can become large enough to account for Mz of HD 149026b. 

In principle, it is possible for a 70M e core to form at ~ 0.05 AU through accretion of many 
Earth-size to Uranian-size cores that migrate close to the parent star. In the vicinity of the parent 
star, gas accretion should be truncated at a relatively small mass (see eq. [14]). We may need 
to carry out N-body simulations of the migrating protoplanets in order to examine whether they 
can pass through resonant trapping (or shepherding) to merge or not. However, we should keep 
in mind that incorporation of type-I migration without any conditions causes severe inconsistency 
with observed data of extrasolar planets and our Solar System. If we rely on type-I migration 
model, we need to clarify the condition for the occurrence of type-I migration at the same time. 

5. SUMMARY AND CONCLUSIONS 

The high density of the close-in extrasolar giant planet HD149026b recently discovered by Sato 
et al. (2005) challenges theories of planet formation. In this paper, we have attempted to derive 



-17- 



robust constraints on the planet's composition and infer possible routes to explain its formation. 

We have first simulated the evolution of HD 149026b more extensively than previous workers 
(Sato et al. 2005; Fortney et al. 2006) and confirmed that the planet contains a substantial amount 
of heavy elements. Preferred values of the total mass of heavy elements are 50-80 M ffi (section 2.1), 
which is consistent with the previous calculations. We showed that the results are unchanged for 
heavy elements located in the central core, or distributed inside the envelope, provided they remain 
deeper than the external radiative zone. In the event of a significant enrichment of the outer 
layers, slightly higher values of heavy elements content are possible (section 2.3). In order to derive 
minimum values of the mass of heavy elements, we have explored the possibility that the planet 
was stored in a relatively cold environment for some time before migrating near to the planet. This 
strict minimum is ~ 35M0, but is regarded as unlikely because it requires a late migration and no 
reheating of the planet by tidal circulization (section 2.2). 

We have then investigated the possibility of subcritical core accretion as envisioned for Uranus 
and Neptune to account for the small envelope mass as well as the large core mass. In principle a 
large core of 50-80 M e can be formed by subcritical core accretion. However we have found it very 
unlikely for at least two reasons: (i) A subcritical core accretion results in a ratio of the core mass to 
the total mass above ~ 0.7 (section 3.1), whereas our evolution calculations showed such a high ratio 
to be possible in a very limited range of parameters (see Table 1); (ii) The subcritical formation of 
a 50-80 M ffi core requires an extremely massive or metal-rich disk with dust surface density 30-50 
times the values obtained for the minimum mass solar nebula (section 3.2). A reasonably massive 
and/or metal-rich disk can form cores of at most ~ 30M ffi far from the parent star. Those facts 
require us to consider (i) the migration of the planet, (ii) the supply of heavy elements to the planet 
during or after the gas accretion phase, and (iii) a limited supply of disk gas or loss of the envelope 
gas to account for the properties of HD149026b. 

In section 4.1, we have discussed how the heavy elements can be delivered to the planet during 
or after the gas accretion phase according to these scenarios. An efficient delivery during the gas 
accretion phase needs to be re-investigated in much more details, because the shepherding tends 
to prevent the planet from accreting planetesimals (section 4.1.1). On the other hand, scattering 
of planetesimals/planets by one or several outer giant planets was shown to lead to an efficient 
accretion by a close-in giant planet, and is a promising explanation for the formation of metal-rich 
planets (section 4.1.2). 

The relatively small amount of hydrogen and helium present in HD149026b is also to be 
explained by the relatively slow formation of the planet in a low-density environment (section 4.2.1). 
Another possibility is related to erosion of the envelope during giant impacts. A reasonable impact 
velocity of 2-2.5 times the two-body surface escape velocity was found to result in a substantial 
loss of the envelope gas, while solid cores are probably merged (section 4.2.2). 

In summary, we can propose at least two scenarios for the origin of HD149026b. 
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• A giant planet with a core of ~ 30M e forms far from its parent star in a relatively massive 
and/or metal-rich disk. Then it moves to the vicinity of the parent star through type-II 
migration (Lin et al. 1996). In the outer regions additional more than three giant planets 
form, which may be likely in massive/metal-rich disks (Ida & Lin 2004a, b). They starts orbit 
crossing and the innermost one is temporally detached from the outer two. The perturbation 
from the outer two secularly pumps up the eccentricity of the innermost one (Nagasawa et 
al. 2006). The associated secular decrease of its pericenter distance results in collision with 
the close-in planet. The impactor also has a core of ~ 30M®. Their solid cores are merged to 
form a core of ~ 6OM0, while envelope gas is severely eroded to MjjHe ~ 5OM0. The orbital 
eccentricity of the merged body is damped by tidal interaction with the parent star. 

• A ~ 30M® core forms in an evaporating disk that was originally massive and/or metal rich. 
The core thus becomes supercritical, but the accretion of the envelope is limited by viscous 
diffusion in an evolved, extended gaseous disk (see Guillot & Hueso 2006). The planet 
subsequently accretes heavy elements through the delivery of planetesimals and/or planets, 
either because of migration due to interactions with the gas disk or to secular perturbations 
with a massive outer planet. 

Above two scenarios are different from each other in that the former leaves at least one outer 
giant planet while the latter does not necessarily need any outer giant planets. Long-term radial 
velocity observations of HD 149026b are thus needed to determine the presence of other planets in 
the system. 

We thank the anonymous referee for fruitful comments. This research was supported by Min- 
istry of Education, Culture, Sports, Science and Technology of Japan, Grant-in-Aid for Scientific 
Research on Priority Areas, "Development of Extra-solar Planetary Science" . 
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Fig. 1. — Evolution of the radius of a 0.36 Mj planet made of a central core of 30 to 75 and a 
gaseous envelope of solar composition. The planet is irradiated so that T e * q = 1740 K (see text). The 
core is assumed to be made of ices (plain lines) or rocks (dashed lines). The results are compared 
to the observational constraints on the age and radius of HD149026b (black box). 
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Fig. 2. — Evolution of the radius of a 0.36 Mj planet. The plain line corresponds to the case of 
a 60 ice core with a solar composition envelope shown in Fig. 1. The dotted line represents 
the result of the calculation with a 30 M® ice core and 30 M ffi of ices in the envelope; opacities 
were calculated for a solar-composition mixture. The dashed line is the same calculation, but with 
an opacity that is artificially increased by 30 to account for the presence of heavy elements in 
the envelope. The results are compared to the observational constraints on the age and radius of 
HD 149026b (black box). 
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Fig. 3. — The critical core mass as a function of core accretion rate for different choices of distance 
from the parent star and local density of disk gas. 
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Fig. 4. — The ratio of core mass to planetary total mass as a function of core mass normalized by 
the critical core mass. The dotted line represents the result for M c of 1 x lO -5 M0/yr at 0.1 AU, 
the dashed one for M c of 1 x lO -5 M0/yr at 1 AU, and the solid one for M c of 1 x 10~ 3 M e /yr at 
5.2 AU. 
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critical core mass (Earth mass) 

Fig. 5. — Timescale for gas accretion (solid line) is shown as a function of the critical core mass 
together with time to the critical core (dashed lines). The timescale for gas accretion is that for an 
e-fold increase in envelope mass just after the critical core mass is attained (see text for the exact 
definition). The circles, triangles, and crosses are for 0.1, 1, and 5.2 AU in the minimum-mass solar 
nebula. The vertical part of the dotted curve with circles corresponds to the horizontal part of the 
corresponding line in Fig. 3. 
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Fig. 6. — The conditions for subcritical core accretion. The subcritical core accretion is possible 
in the regions over the solid, dotted, and dashed lines for 30, 45 and 70M ffi cores, respectively. 
Corresponding to HD149026, M* = 1.3M and [Fe/H]=0.36 (/ d // 6 = 2.3) are assumed. Note 
that the discontinuities of the curves come from sudden increases in disk surface density of dust 
component at the snow line. 
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Fig. 7. — Probability of outcomes of encounters between the inner giant planet of 0.5Mj in a circular 
orbit at 0.05 AU and a planetesimal (test particle) (panel a) or another giant planet of 0.5Mj (panel 
b) in a nearly parabolic orbit with initial semimajor axis a = 1 AU. For each initial q, 100 cases 
with random angular distributions are calculated until 100 Keplerian periods at 1 AU. Collision 
probability with the inner planet (Pcoi)j that with the parent star (-P*i), and ejection probability 
(P e jc) are expressed by solid, dotted, and broken lines. Here, the physical size of planets are 
determined with internal density 1 gem -3 ; that of the parent star is set as 0.01 AU. 
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Fig. 8. — The loss fraction of mass after a head-on collision between two Jupiter-mass planets 
composed of ideal gas with the adiabatic exponent of 2. The loss fraction of mass initially located 
inside 1/2 (dashed line) and 1/4 (dotted line) of the initial planetary radius are plotted to compare 
with total loss fraction (solid line). The loss fraction is defined as the mass fraction of the SPH 
particles whose velocity exceeds the local escape velocity and distance from the center of mass 
frame exceeds 5 times of the initial planetary radius. 



